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The  NIST  Algorithm  Testing  System  (or  ATS)  is  a software  package  for  testing  the 
performance  of  geometric  fitting  software  used  in  metrology  systems'.  One  of  the  functions  of 
the  ATS  is  to  generate  test  data  simulating  point  measurements  taken  from  manufactured  parts. 
Data  sets  are  generated  by  evaluating,  according  to  user-specified  sampling  plans,  models  of  parts 
with  simulated  form  errors.  The  ATS  also  provides  a simple  model  of  point  measurement  errors. 
The  purpose  of  this  report  is  to  document  the  form  error  models  used  in  the  ATS. 

The  ATS  uses  a parametric  representation  of  geometry.  A form  error  is  simulated  by 
computing,  from  the  parametric  coordinates  of  a point  on  an  ideal  geometric  element,  the 
cartesian  coordinates  of  a form  error  vector  at  that  point.  This  error  vector  is  added  to  the  ideal 
point  vector  to  generate  a point  on  a simulated  surface  with  form  error.  A simulated 
measurement  point  is  finally  generated  by  adding  a random  measurement  error  vector.  In 
summary,  a point  generated  by  the  ATS  data  generator  is  of  the  form: 

p(w,v)  = pfu,v)  + App(M,v)  + Apj^ 

where  (m,v)  are  the  parametric  coordinates  of  the  point  (according  to  the  sampling  plan),  Po(m,v) 
is  the  point  on  the  ideal  geometry,  App(M,v)  is  the  form  error,  and  Ap^^  is  the  measurement  error 
at  that  point.  (The  above  description  applies  to  surface  geometries.  Curve  geometries  are 
similar,  but  only  use  a single  evaluation  parameter.) 

The  ATS  simulates  different  kinds  of  form  errors  for  different  ideal  geometries.  However, 
a random  form  error  is  supported  for  all  curve  and  surface  types.  This  consists  of  a displacement 
perpendicular  to  the  ideal  geometry  by  some  random  distance.  Random  form  errors  are  not 
reproducible:  evaluating  the  same  surface  at  the  same  parametric  coordinates  will  produce  a 
different  point.  For  other  than  random  form  errors,  the  ideal  geometries  and  corresponding  form 
error  types  are  as  follows: 


Geometry 

Form  Error 

Description 

Line 

Sine 

A sinusoidal  oscillation  perpendicular  to  a line. 

The  oscillation  takes  place  in  a plane  determined 
by  the  ideal  line  definition. 

Step 

A discontinuity  in  a line.  After  the  step,  the  line 
continues  in  the  same  direction,  but  displaced  from 
the  original  line.  The  displacement  is  in  a plane 
determined  by  the  ideal  line  definition. 

Bend 

A change  in  direction  of  a line.  At  the  bend,  the 
line  is  continuous,  but  the  direction  of  the  line  is 
discontinuous  The  new  line  lies  in  a plane  deter- 
mined by  the  ideal  line  definition. 

An  overview  report  on  the  ATS  is  in  preparation. 


Geometry 

Form  Error 

Description 

Circle 

Sine 

A sinusoidal  oscillation  perpendicular  to  a circle,  in 
the  plane  of  the  circle.  To  keep  the  curve  simple, 
the  amplitude  of  the  oscillation  should  be  less  than 
the  radius  of  the  circle. 

Step 

A discontinuity  in  the  radius  of  the  circle  between 
two  polar  angles. 

Plane 

Sine 

A two-dimensional  sinusoidal  oscillation  about  a 
plane.  The  oscillation  is  perpendicular  to  the 
plane. 

Step 

A discontinuity  in  a plane  at  a line  in  the  plane. 

One  one  side  of  the  line  of  discontinuity,  the  plane 
is  displaced,  but  parallel  to,  the  plane  on  the  other 
side. 

Bend 

A sharp  bend  or  crease  in  a plane  along  a line.  On 
both  sides  of  the  line,  the  plane  has  local  ortho- 
gonal coordinate  systems.  The  normals  are  differ- 
ent. 

Sphere 

Sine 

A two-dimensional  sinusoidal  oscillation  perpendic- 
ular to  a sphere.  To  keep  the  surface  simple,  the 
maximum  amplitude  of  the  oscillation  should  be 
less  than  the  radius  of  the  sphere. 

Cylinder  or 
Cone 

Axial  sine 

A sinusoidal  oscillation  of  the  axis.  The  oscillation 
takes  place  in  a plane  containing  the  ideal  axis. 

(Refer  to  the  sine  error  for  a line.) 

Radial  sine 

A sinusoidal  oscillation  of  the  radius  at  each  cross 
section  perpendicular  to  the  axis.  (Refer  to  the 
sine  error  for  a circle.) 

Axial  radius  sine 

A sinusoidal  oscillation  of  the  radius  along  the 
axis. 

Axial  step 

A discontinuity  of  the  axis.  (Refer  to  the  step  error 
for  a line.) 

Radial  step 

A discontinuity  in  the  radius  at  each  cross  section 
perpendicular  to  the  axis.  (Refer  to  the  step  error 
for  a circle.) 

Axial  radius  step 

A discontinuity  in  the  radius  along  the  axis. 

The  tables  on  the  following  pages  summarize  the  form  error  models  used  in  the  ATS. 
There  is  one  table  for  each  type  of  ideal  geometry  supported  by  the  ATS.  Each  table  presents 
the  types  of  form  errors  supported  for  that  geometry,  the  parameters  that  define  each  error  type, 
a name  for  each  form  error  parameter,  and  an  equation  for  the  error  vector  as  a function  of 
parametric  coordinates  of  points  on  the  ideal  surface. 
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Throughout  these  tables,  the  following  special  functions  are  used: 


• The  unit  step  function; 


jc<0 


x>0 


• The  positive  ramp  function; 


f 0 x<0 


TU) 


X jc>0 


• The  function  U(a,B):  a uniformly-distributed  random  variate  in  the  range  [a,B]. 

The  tables  also  assume  that  each  geometry  has  a local  cartesian  coordinate  system,  which 
is  used  to  evaluate  each  form  error  vector.  The  local  coordinate  system  is  described  briefly  at 
the  end  of  each  table.  The  geometry  of  the  ATS  is  based  on  STEP^.  More  complete 
information  about  geometry  representations  used  in  the  ATS  is  in  preparation. 


^ STEP  is  an  acronym  for  an  emerging  international  standard,  the  Standard  for  the  Exchange  of  Product  Model  Data. 
The  geometry  of  the  ATS  is  based  on  Part  42  of  STEP,  specifically,  the  following  document  of  ISO/TC  184/SC  4: 
Product  Data  Representation  and  Exchange — Part  42:  Integrated  Resources:  Geometric  and  Topological  Representation, 
ISO  Committee  Draft  10303-42,  June  21,  1991. 
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LINES 


Error 

Parameter 

Parameter  name 

Equation 

Sine 

A 

amplitude 

Asm(Ju+p)n^ 

/ 

frequency 

P 

phase 

Step 

A 

amplitude 

ASiu-u^n^ 

Mo 

location 

Bend 

n, 

new  direction 

T(M-Mo)(n,-n) 

Mo 

location 

Random 

A 

amplitude 

U(-A,A)n^ 

The  local  coordinate  system  consists  of  n,  the  direction  of  the  line. 


In  the  above  equations. 


n-z  I > 0.9 
n*z  I < 0.9 

is  the  unit  direction  perpendicular  to  n closest  to  the  +z  axis  (or  +x  axis  if  n is  close  to  z).  Also, 
if  the  line  is  parameterized  from  to  u^. 


W IX  ■ 

max  nun 

is  the  normalized  coordinate  of  the  point  in  the  line  segment  being  modeled.  (In  other  words, 
the  frequency  / indicates  the  total  number  of  wavelengths  of  form  error  oscillation  in  the  line 
segment.) 


nxx 


nxz 
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CIRCLES 


Error 

Parameter 

Parameter  name 

Equation 

Sine 

A 

amplitude 

Asin(Ju+p)n{u) 

/ 

frequency 

P 

phase 

Step 

A 

amplitude 

A(S(m-Mq)  -S(M-Mj))n(M) 

Mo 

start  angle 

Ml 

end  angle 

Random 

A 

amplitude 

U(-A,A)n(M) 

The  local  coordinate  system  consists  of  and  n2,  two  perpendicular  directions  in  the  plane  of 
the  circle. 

In  the  table, 

n(w)  = (cosi/)nj +(sinM)n2 

is  the  unit  direction,  away  from  the  center,  perpendicular  to  the  circle  at  u. 
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PLANES 


Error 

Parameter 

Parameter  name 

Equation 

Sine 

amplitude  in  u 

[Asinifu-^pJ 

+^/in(/v+p^)]n3 

fu 

frequency  in  u 

Pu 

phase  in  u 

Av 

amplitude  in  v 

/v 

frequency  in  v 

p. 

phase  in  v 

Step 

A 

amplitude 

AS(Msina  -vcosa-uOn^ 

d 

distance 

a 

orientation 

Bend 

/ 

"3 

new  normal 

T(«sina  -vcosa  -/)  • 

1 / 1 ("3  ^"3) 



|n3^"3  1 

d 

distance 

a 

orientation 

Random 

A 

amplitude 

U(-A,A)n3 

The  local  coordinate  system  consists  of  nj,  02,  and  n3,  three  mutually  perpendicular  unit  vectors, 
with  n3  perpendicular  to  the  plane. 


Also,  if  the  plane  is  parameterized  from  to  in  u,  and  similarly  for  v. 


u-u 

a = 27t  - . , 


u -u 

max  min 


V-V  . 

V = 27t — 

V -V 
max  min 

are  the  normalized  coordinates  of  the  point  within  the  planar  patch  being  modeled.  (In  other 
words,  the  frequencies /„  and/,  indicate  the  total  number  of  wavelengths  of  form  error  oscillation 
in  the  u and  v extents.) 
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SPHERES 


Error 

Parameter 

Parameter  name 

Equation 

Sine 

amplitude  in  u 

[Asinifu^pJ 

+A^smifv+pJ]n(u,v) 

fu 

frequency  in  u 

Pu 

phase  in  u 

amplitude  in  v 

/v 

frequency  in  v 

p. 

phase  in  v 

Random 

A 

amplitude 

U(-A,A)n(M,v) 

The  local  coordinate  system  consists  of  nj,  n2,  and  Oj,  three  mutually  perpendicular  unit  vectors. 
In  the  above  equations: 

n(w,v)  = (cosw)(sinv)nj +(sinM)(sinv)n2 +(cosv)n3 
is  the  unit  direction,  away  from  the  center,  perpendicular  to  the  sphere  at  (m,v). 
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CYLINDERS  and  CONES 


Error 

Parameter 

Parameter  name 

Equation 

A 

amplitude 

Axial  Sine 

/ 

frequency 

AsinC/v+p)n(Mo) 

P 

phase 

Mo 

reference  angle 

A 

amplitude 

A sm(fu+p)n(u) 

Radial  Sine 

/ 

frequency 

P 

phase 

A 

amplitude 

Asin(/v+p)n(M) 

Axial  Radius 
Sine 

f 

frequency 

P 

phase 

A 

amplitude 

AS(v-Vo)n(Mo) 

Axial  Step 

Vo 

location 

Mo 

reference  angle 

A 

amplitude 

A(S(m-Mq)  -S(M-Mj))n(M) 

Radial  Step 

Mo 

start  angle 

Ml 

end  angle 

Axial  Radius 

A 

amplitude 

AS(v-Vo)n(M) 

Step 

Vo 

location 

Random 

A 

amplitude 

U(-A,A)n(M) 

The  local  coordinate  system  consists  of  n,,  02,  and  n3,  three  mutually  perpendicular  unit  vectors, 
with  n3  directed  along  the  axis  of  the  cylinder  or  cone. 

In  the  above  equations, 

n(M)  = (cosw)nj +(sinw)n2 

is  the  unit  direction,  away  from  the  axis,  perpendicular  to  the  axis  at  u (for  all  v).  Note  that,  in 
the  case  of  a cone,  the  form  error  is  not  perpendicular  to  the  surface. 


Also,  if  the  cone  is  parameterized  from  to  in  v (axially), 

v-v 

V = 2k — 

V -V 

max  rmn 

is  the  normalized  v coordinate  of  the  point  within  the  axial  segment  being  modeled.  (In  other 
words,  the  frequency  / indicates  the  total  number  of  wavelengths  of  form  error  oscillation  along 
the  axis.) 
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